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ABSTRACT 

The pulse profiles of the accreting X-ray millisecond pulsar SAX J1808.4— 3658 at different 
energies are studied. The two main emission component, the black body and the Comptonized 
tail that are clearly identified in the time-averaged spectrum, show strong variability with 
the first component lagging the second one. The observed variability can be explained if the 
emission is produced by Comptonization in a hot slab (radiative shock) of Thomson optical 
depth ^ 0.3-1 at the neutron star surface. The emission patterns of the black body and the 
Comptonized radiation are different: a "knife"- and a "fan"-like, respectively. We construct 
a detailed model of the X-ray production accounting for the Doppler boosting, relativistic 
aberration and gravitational light bending in the Schwarzschild spacetime. We present also 
accurate analytical formulae for computations of the light curves from rapidly rotating neutron 
stars using formalism recently developed by Beloborodov (2002). Our model reproduces well 
the pulse profiles at different energies simultaneously, corresponding phase lags, as well as 
the time-averaged spectrum. We constrain the compact star mass to be bounded between 1.2 
and 1.6M0. By fitting the observed profiles, we determine the radius of the compact object 
to be i? - 11 km if M = I.6M0, while for M = 1.2Mo the best-fitting radius is 6.5 km, 
indicating that the compact object in SAX J 1808.4—3658 can be a strange star. We obtain a 
lower limit on the inclination of the system of 65°. 

Key words: accretion, accretion discs - methods: data analysis - pulsars: individual (SAX 
J1808.4-3658)- X-rays: binaries 



1 INTRODUCTION 

Discovery of coherent oscillations with frequencies in the 300-600 
Hz range during X-ray bursts from a number of low-mass X-ray bi- 
naries (see Strohmayer & Bildsten 2003 for a review) has triggered 
the efforts to use the information on the amplitude of variability and 
the folded pulse shape to put constraints on the compactness of the 
neutron star and thus its equation of state as well as the emission 
pattern from the neutron star surface (e.g. Miller & Lamb 1998). 
However, a rather limited photon statistics does not allow to reach a 
high accuracy (Strohmayer et al. 1997; Nath, Strohmayer & Swank 
2002; Muno, Ozel & Chakrabarty 2002). 

The millisecond coherent pulsations discovered in the per- 
sistent emission of the four sources: SAX J1808.4— 3658 (with 
the period P — 2.5 ms; Wijnands & van der Klis 1998a), 
XTE J175 1-305 (P = 2.3 ms; Markwardt et al. 2002), XTE 
J0929-314 (P = 5.4 ms; Galloway et al. 2002) and XTE 
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J1807-294 (P = 5.25 ms; Markwardt, Smith & Swank 2003) 
allow one to increase the statistics by folding the profile over a 
longer observational period (days rather than seconds in the X-ray 
burst oscillations). Thus, for example, SAX J1808.4— 3658 showed 
coherent pulsations with rms amplitude of 5-7 per cent during its 
April 1998 outburst (Cui et al. 1998) with almost constant shape of 
its energy spectrum (Gilfanov et al. 1998) and very similar pulse 
profiles. The pulse profiles of SAX J1808.4— 3658 show strong en- 
ergy dependence with soft photons lagging behind the hard ones 
(soft phase lags, see Cui et al. 1998). By analysing the phase- 
resolved spectra, Gierlinski, Done & Barret (2002) showed that this 
results from the fact that the two main spectral components, a soft 
black body and a hard Comptonized, vary out of phase with the first 
lagging the last one. While the variations of the black body flux can 
be described by a single sine wave, the pulse of the hard Comp- 
tonized component is strongly skewed. Gierlinski et al. (2002) sug- 
gested that the Doppler boosting can play a role in changing the 
shape of the profile. 

The radiation pattern of standard accreting X-ray pulsars is in- 
fluenced by strong magnetic field B ~ 10^^ G. On the other hand, 
the magnetic field of the accreting millisecond pulsars in low mass 
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X-ray binaries is much weaker, B = 10*-10^ G (see e.g. Wijnands 
& van der Klis 1998a; Psaltis & Chakrabarty 1999), and does not 
affect significantly properties of the emitted radiation. Thus these 
sources can also serve as laboratories for studying radiative pro- 
cesses at the surface of weakly-magnetized neutron stars. 

In this paper, we construct a detailed model for the X-ray 
emission from the surface of a rapidly rotating neutron star account- 
ing for relativistic effects. We compare the model with the data on 
SAX J1808.4— 3658 and put constraints on the inclination of the 
system, the position of the emitting region relative to the rotational 
pole and its emission pattern as well as the stellar radius. The data 
used for the analysis are described in §|2| The model and useful an- 
alytical formulae for light curve calculations are presented in § |3| 
The main results are given in §|4| The discussion and the summary 
are given in §|5|and §|6| respectively. 




Figure 1. The geometry of the problem. Dotted curve shows the photon 
trajectory. 



2 DATA 

We study the data obtained by the Rossi X-ray Timing Explorer 
(RXTE) during the April 1998 outburst of SAX J1808.4-3658, 
using FTOOLS 5.2. To improve the statistics we average the data 
between 1998 April 11 and 29 (resulting in 118 ks of data). For 
spectral fitting we extract the Proportional Counter Array (PCA) 
spectra from all five detectors, top layer only, and use the data in 
3-20 keV energy band. We also extract spectra from High-Energy 
X-ray Timing Experiment (HEXTE) in the 20-150 keV band, from 
clusters and 1. We constructed the energy-dependent pulse pro- 
files from the PCA data (all detectors, all layers) by following the 
procedure described in Cui et al. (1998) and Gierlinski et al. (2002) 
correcting the photon arrival times for orbital motions of the pulsar 
and the spacecraft. For pulse profile fitting we created two folded 
light curves in 16 phase bins, for energies 3^ keV and 12-18 keV. 
The statistical uncertainties are about the same as uncertainties in 
the background, resulting in errors of 0.2 per cent and 0.3 per cent, 
respectively, in the count rate in the two considered energy bands. 



3 MODEL 

3.1 Calculational method 

Accreting matter following the magnetic field lines close to the 
neutron star is stopped in the very vicinity of the surface by ra- 
diation producing radiation dominated shock (Basko & Sunyaev 
1976; Lyubarskii & Sunyaev 1982). For the source luminosity of 
a few percent of the Eddington luminosity (Gilfanov et al. 1998; 
Gierliriski et al. 2002), the vertical (i.e. along the radial direction) 
extension of the shocked plasma is smaller than a characteristic hor- 
izontal size and certainly smaller than the stellar radius. The hard 
X-rays produced in the shock can irradiate the surrounding stellar 
surface, so that the black body emission region can cover a some- 
what larger area. For our calculations we assume that all photons 
originate at the stellar surface (i.e. the height of the emitting region 
is set to zero). 

We assume a circular spot and consider two extreme cases 
for the relative positions of the black body and the Comptonizing 
regions: (1) a homogeneous slab, i.e. the hot Comptonizing layer 
covering the whole black body spot, (2) a point-like Comptonizing 
region in the centre of a black body region. 

In order to compute the light curve as observed by a distant 
observer, we first specify the radiation spectrum and angular depen- 
dence in the frame co-rotating with the star. We then make Lorentz 



transformation to obtain the radiation intensity in the non-rotating 
frame close to the stellar surface. At the final step, we follow pho- 
ton trajectories in the Schwarzschild spacetime to infinity. 

Let us now consider a star with an azimuthally symmetric sur- 
face radiation intensity l'{E',a') that can vary over the surface, 
where a' is the angle between an emitted photon and the local nor- 
mal to the stellar surface. Let dS' be a surface element (spot) at 
colatitude 9 (see Fig. Q for the geometry). The primed quantities 
correspond to the frame co-rotating with the spot. Let fe and n be 
unit vectors pointing from the star centre towards the observer and 
the spot, respectively, and i be the inclination of the spin axis. The 
inclination of the spot varies periodically 

cost/; — k ■ n — cos i cos + sini sin cos tj), (1) 

where the phase — 2111/1, with u being the pulsar frequency and 
t = is chosen when the spot is closest to the observer. 

We compute the original direction of the photon feo near the 
stellar surface (which is transformed to k at large distance from the 
star) assuming Schwarzschild geometry where photon orbits are 
planar: 

fco = [sin a fe + sin(?/' — a) n]/ sinip, (2) 

where a is the angle between fco and n, i.e. cos a = ko ■ n. 

The relation between a and tp (i.e. light bending) can be ob- 
tained by standard techniques (Pechenick, Ftaclas & Cohen 1983): 

J r'^ \p- V r I ' 

where the impact parameter 

6 — — ^ - sin Q, (4) 

Tg = 2GM/(? is the Schwarzschild radius, M is the mass and R 
is the radius of the compact star. 

The observed flux at energy E is AFe ~ I{E)d^l, where 
I{E) is the radiation intensity at the infinity and dQ is the solid 
angle occupied by dS' on the observer's sky. The solid angle can 
be expressed through the impact parameter 

dfl^bdbdip/D^, (5) 

where D is the distance to the source and ip is the azimuthal an- 
gle corresponding to rotation around fc. The impact parameter b 
depends on t/j only, but not on ip. 

The apparent area of the spot as measured by photon beams 
in the non-rotating frame near the stellar surface is dS = 5 dS' 
(see Terrell 1959; Lightman et al. 1975; Ghisellini 1999) and the 
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Figure 2. Angular distribution of the radiation flux (normalized as J fil{ij,)dfi = 1) escaping from the top of an electron scattering slab of Thomson optical 
depth Tes = 0.3 (left) and Tcb = 0.7 (right). Incident radiation from the bottom is a black body (intensity independent of the zenith angle arccos fi). Com- 
putations follow procedure described in Sunyaev & Titarchuk (1985). Different scattering orders are shown and marked by numbers. In the Comptonization 
process, scatterings also shift photons along the energy axis so that at a given energy one scattering order dominates (see Sunyaev & Titarchuk 1985; Poutanen 
& Svensson 1996). 



relation between a and a' is described by the relativistic aberration 
formula (for motions parallel to the spot surface) 



cos a' = 5 cos a, 



(6) 



where S is the Doppler factor. Thus the spot area projected on to 
the plane perpendicular to the photon propagation direction, i.e. a 
photon beam cross-section, is Lorentz invariant (see e.g. Lightman 
et al. 1975; Lind & Blandford 1985): 



dS cos a — dS' cosq'. 



(V) 



Representing dS — d cos dip and using equations ^4) and Q 
we get from jsj 



dn 



dS' cos a' 



dcosa 



(8) 



l — Tg/R dcos?/) 
The Doppler factor can be expressed as follows 

5 = l/[7(l-/3cose)], (9) 

where 7 = (1 — 0^)^^^'^ and 13 — v/c is the spot velocity as 
measured in the non-rotating frame, 

2tvR V 



/3 = 



c ^\-r,IR 



sin d — /3cq sin ( 



(10) 



Here /3cq is the velocity at the equator and ^ is the angle between 
the spot velocity and fco. The pulsar frequency is corrected here for 
the redshift \/l — ^a/R- Using equation J2j it is easy to show that 



/3 , sin a /3 , 
cosC = -J • fco = -. — 7-j • fe 

p simp p 



sin 7/; 



(11) 



In the case of a weak field and negligible bending the solid 
angle occupied by the spot is dfl = dS' cos a' /D^. In reality, this 
formula can be applied even when bending is significant since the 
relation between cos a and cos is close to linear (Beloborodov 
2002), 



rg/R + (1 — r-g/R) cosi/', 



(12) 



for a star with R > 2r\j, so that d cos a/d cos 1/) ~ 1 — rg/7?. 
The radiation intensity observed at the infinity is I{E) 



(1 - rg/Rf^^Io{EQ, a), where E/Eq = ^1 - rJR is the red- 
shift (Misner, Thorn & Wheeler 1973) and the energy _Z?o is mea- 
sured in the non-rotating frame near the stellar surface. (In the ac- 
tual calculations, one can neglect the redshift factor in all formulae 
for photon energies, since it is the same for the photons originat- 
ing at the stellar surface.) Now we can make Lorentz transforma- 
tion to the spot rest frame where we know the angular and energy 
distribution of the radiation field I'{E' , a'). The intensities are re- 
lated via Jo {Eo , a) = S'-^I' {E' , a) with the Doppler shift given by 
Eo = 5 E'. 

Combining equations above we get the expression for the ob- 
served flux: 



dFE 



-fyi ^ ^ [E ,a) 



d cos a 
dcos?/; 



dS' cos a', 



(13) 



with the visibility condition cos a' > 0. We account thus here for 
the special relativistic effects (Doppler boosting, relativistic aberra- 
tion), the gravitational redshift and light bending in Schwarzschild 
geometry. In order to evaluate the flux at a given phase 0, one first 
computes ip using equation Q, then inverting js} gets a and sub- 
stitutes the results into equations <1 1> and j9j to obtain the Doppler 
factor which is used later to compute a' from equation j6|. The 
light curve from the antipodal spot can be easily obtained by re- 
placing in all formulae 9 ^ n — 8 and <j} ^ n + (j). We note 
here that the formalism outlined above neglects the time delays due 
to different photon paths (see Pechenick et al. 1983) which can be 
easily incorporated. However, they are negligible comparing to the 
pulsar period above a few ms. 

For a spot of the finite size, the observed flux can be obtained 
by integrating over the (visible) spot area. On the other hand, if we 
know from the observations the black body flux _Fbb, its tempera- 
ture Too and the distance to the source D, one can estimate the size 
of the (black body) emitting region. The "observed" radius of the 
spot rspot.oo = {FhhD'^ /o'ssT^y^'^ should be corrected for the 
effects of gravitation and orientation (using a procedure described 
in 13.31 see also Psaltis, Ozel & DeDeo 2000) to obtain the actual 
spot area. 

In our calculations of the model light curves for SAX 
J1808.4-3658 we use the distance D = 2.5 kpc (in't Zand et 
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Figure 3. Time-averaged spectrum of SAX J1808.4— 3658 as observed 
by RXTE. The model spectrum (solid curve) consists of a black body 
(dashes/dots), thermal Comptonization (dots), and Compton reflection 
with the iron line (dashes). Thermal Comptonization component (model 
THCOMP in XSPEC; Zdziarski, Johnson & Magdziarz 1996) has a best- 
fitting photon spectral slope V = 1.86 (the electron temperature was fixed 
at 90 keV) and the Compton reflection amplitude C/(27r) = 0.09. The 
black body modeled by BBODYRAD from XSPEC has the best-fitting tem- 
perature fcTbb = 0.66 keV and the radius of the circular emitting region 
''spot.oo = 2.4 Ian (for the source distance D = 2.5 kpc). 

al. 2001), Too = 0.66 keV, and rspot.oo = 2.4 km (see Fig.|3l 
and Gierlinski et al. 2002). Since the observed profiles from SAX 
J1808.4— 3658 have almost sinusoidal shapes, they cannot be pro- 
duced by two antipodal spots. Therefore, we compute the light 
curves from the primary spot (closest to the observer) only. We 
also neglect multiple images since neutron or strange stars have 
sizes larger than l.Srg (see e.g. Lattimer & Prakash 2001). 

3.2 Incident spectrum and its angular distribution 

In order to compute a light curve Fe {<j>) at a given energy E, one 
needs to specify the intrinsic radiation pattern in the spot rest frame 
I'{E',a'). In a slab geometry of the emitting region, the black 
body photons come from the bottom of the hot Comptonizing slab. 
The black body intensity transmitted through the slab of Thomson 
optical thickness r = n^aTh is oc exp(— r//i), where 

arccos is the zenith (polar) angle measured from the slab nor- 
mal. Tie is the electron concentration, ax is Thomson cross-section, 
and h is the height of the slab. The radiation flux from a unit area 
is then strongly peaked along the normal direction (a "knife"-like 
emission pattern, see curves marked by at Fig.|2j- The scattered 
radiation is expected to have a very different angular distribution, 
which depends on the scattering order (and thus on photon energy). 
The radiation flux does not peak anymore in the normal direction, 
but instead peaks at some intermediate angle (a "fan"-like emis- 
sion pattern, see solid curves marked by oo in Fig.|2j- It approaches 
the asymptotic distribution 7oo (t, /x) in a few scatterings (the exact 
number is a function of the optical depth). We compute the angular 
distribution of radiation following procedure described in Sunyaev 
&Titarchuk(1985). 

The intensity of the black body radiation in the frame co- 
moving with the spot can be represented in the following way: 

/^^(a') oc exp(-r/ cos a')/bb(-B'\/l^7^/E), (14) 

where fhh{E) is the observed time-averaged black body flux (i.e. 
the best-fitting BBODYRAD model spectrum, see Fig.|3j. We as- 



sume that the energy and angular dependences of the high energy 
(scattered many times) radiation can be separated as 

7|?,(a') (X I^{T,a')fsciE'^l-rJR). (15) 

Here fsc{E) is the observed, time-averaged energy dependence of 
the Comptonized flux (i.e. the best-fitting THCOMP model spec- 
trum, see Fig.|3}. 

Since the main spectral components of the time-averaged 
spectrum of SAX J1808.4— 3658 are a power-law-like Comp- 
tonized component (with weak Compton reflection) extending at 
least to 200 keV and a black body contributing about 30 per cent to 
the flux in the 3-5 keV region, (see Fig.|3|and Gilfanov et al. 1998; 
Gierliiiski et al. 2002), we compute the total observed flux as a 
function of phase Fe{4>) at every energy E summing up the contri- 
butions from the black body and the Comptonized components. We 
compute them independently and renormalise so that their phase- 
averaged values are equal to the observed fluxes /bb(S), /sc(-E) 
(see Fig.|3}. 

If the black body spot is not covered completely by a Comp- 
tonizing layer, some black body photons do not enter the Comp- 
tonizing region and escape to the observer directly. In that case, the 
"effective" r that describes the mean (i.e. averaged over the spot) 
angular dependence of the black body radiation can be smaller 
than the actual optical depth in the Comptonizing region Tcs that 
describes the angular dependence of the Comptonized radiation. 
Since the exact geometry is unknown, we use r and Tcs as two inde- 
pendent parameters for the light curve computations. The emission 
pattern for the black body radiation is given by equation <14t while 
the Comptonized radiation is described by equation <15> where in- 
stead of r we use Tcs as a parameter for the data fitting. These 
patterns are referred to as model 1 in §|4|and TableQ 

In order to check the robustness of the results and their depen- 
dence on the assumed angular distribution of radiation, we consider 
a second model where instead of loo (tcs , m) we use a linear func- 
tion 1 + a/i with a as a free parameter (limited by —1 from below 
to ensure positiveness of the function). The corresponding fits are 
referred below as model 2. 

Thus, the model parameters are the neutron star mass A/, 
its radius 7?, rotational frequency i/ (= 401 Hz for SAX 
J1808.4— 3658), inclination i, colatitude of the spot centre 6, the 
optical depth r and a parameter determining the angular distribu- 
tion of the Comptonized radiation Tcs or a (plus a free phase factor). 

3.3 Analytical light curves 

Using the formalism developed by Beloborodov (2002) we derive 
here simple analytical formulae for the light curves and the oscil- 
lation amplitudes that could be useful for understanding the main 
physical effects. Let us first consider a small spot and make approx- 
imation to the light bending formulae. The bolometric flux can be 
obtained from equation <13> by integrating over E and using the 
cosine relation <12t which describes well the light bending near a 
star with R > 2rg : 

AF = [1 - rj RfS^l' {a')5 cos a AS' /D'^ (16) 

- ('-S'*'B+('-%)-*]'V)f- 

with the visibility condition cos ijj > /i„ = —rg/{R — rg). (For 
the antipodal spot, one should substitute cos i/j ^ — cos and 
— » — (ji when computing the Doppler factor.) 

It is clear from equation <16> that the (bolometric) black body 
flux observed from a spot at a rapidly spinning star is a factor 
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Figure 4. The bolometric black body and Comptonized fluxes, and the 
Doppler factor 5 as functions of phase. Dashed curves correspond to a 
slowly rotating neutron star. Solid curves give the profiles modified by 
the Doppler boosting and aberration for a neutron star of rotational fre- 
quency 401 Hz. (Other pai'ameters are given in Fig.|5]) Light bending in the 
Schwarzschild geometry is accounted for. Dotted curves (coinciding com- 
pletely with the solid one at the upper panel) show the results of multi- 
plication of the dashed curves by 5^ and 5* for the black body and the 
Comptonized emission, respectively (see eqs 1161 and 1171 ). 



times the flux from a slowly rotating star (see Fig.|4p). Two powers 
of 5 come from the solid angle transformation, one from the en- 
ergy, one from the arrival time contraction, and the fifth from the 
change in the projected area due to aberration. The scattered radi- 
ation, on the other hand, has an angular distribution which differs 
from that of the black body. In the optically thin case, it is closer 
to I' [a') cos a' = /' « const (see Fig.|2j. Equation <16t for the 
observed flux can be rewritten then as 

di^sc = (1 - rJRf5*fAS'/D\ (17) 

One sees that the flux is modified by approximately factor 5* (see 
Fig.|4j)). The angular distribution of the hard radiation (defined in 
our model by parameters Tcs or a) is reflected in the shape of the 
light curve at high photon energies, where Comptonizing emission 
dominates. Since the Comptonization spectrum is close to a power- 
law I{E) oc with the same approximations as above, 
one can easily get the monochromatic flux at a given energy: 

dFlf = (1 - rg/Rf'^5^+''f'{Eo)AS'/D^, (18) 

which varies in a way similar to the bolometric flux if F ~ 2. 

Neglecting the Doppler effect and using equation <16> one can 
relate the variability amplitude of the black body flux to the angles 
i and 6 and the neutron star radius (Beloborodov 2002): 

A = (F„,ax - Fn,in)/(Fn,ax + Fn,in) = U/Q, (19) 

where we defined 

U — (1 — rg/i?) sin z sin 

Q = rg/R+{l-rg/R)cosi cose. (20) 

This peak-to-peak amplitude is \/2 larger than sometimes quoted 
rms amplitude if the light curve is a pure sine wave. 



For a circular spot of an angular radius p emitting as a black 
body, the observed flux can be obtained from equation <16> : 



Fo 



"^j p[Q+^ tan - + Ucos( 



,(21) 



if the spot is always visible, i.e. cos(i + 9 + p) > /ii,. Here the 
intensity is related to the local temperature Jo = usbTq /n. We 
can easily now obtain the oscillation amplitude: 

U 



Q + 



■ tan 



2 P ■ 



(22) 



R 2 

Equation <21> allows us to estimate the size of the spot Rp from the 
equation: 



'■spot,oo = sin^ p 



R 2 



(23) 



if the distance D, the temperature = To{l — rg/R^^^ 
and the observed phase-averaged black body flux Fp = 
cseT'o^J'spot.txj/^^ (and thus rspot.oo) are known from observa- 
tions. For small p, the spot radius is Rp = r'spot.ooQ"^''^. 



4 RESULTS 

We fit the observed pulse profiles in the energy bands 3-4 keV and 
12-18 keV simultaneously. The second band is completely dom- 
inated by the Comptonized radiation and thus gives direct infor- 
mation about its angular distribution. In the lower energy band, 
the black body affects the observed pulse profile. Since the mass 
of the compact star in SAX 11808.4—3658 is not known, we con- 
sider masses between A4 — l.OM© and 1.6Mq. Masses in the 
range of 1.4-1.6Mq are expected for neutron stars in low mass 
X-ray binaries that have accreted material from their binary com- 
panions, lower masses are possible, for example, for strange stars. 
We restrict the inclination to i < 82° to be consistent with the 
absence of the X-ray eclipses (Bildsten & Chakrabarty 2001). The 
best-fitting parameters are presented in TableQ The two sets of er- 
rors correspond to the 90 and 99 per cent confidence intervals with 
Ax^ = 2.71 and 6.63, respectively. 

We see from Table Q that the results obtained with the two 
models 1 and 2 are very consistent with each other. For both mod- 
els, a large inclination is preferred and the lower limit is i > 65° 
(at 99 per cent confidence level) for M > 1.2Mq . The fits become 
much worse at masses higher than 1.6Mq. The estimated stellar 
radius (even in Vg units) grows with the mass and varies between 
~ l-7rg and ~ 2.3rg (~ 2.8rg for model 1). Using simple ana- 
lytical estimations we explain the physical reasons behind this be- 
haviour in ^ \5A\ The radii are consistent with some of the neutron 
star equations of state if M > 1.4M0 . For a smaller assumed mass 
of 1.2M0, the radius is more consistent with the equations of state 
for strange stars, while for even smaller mass of IMq, the obtained 
stellar radius of J? ~ 5 km is too small even for a strange star (see 
ij |5.2l and Fig.Qfor details). 

The colatitude 6 of the spot centre varies approximately as 
oc i?-^ Alar ger 6 is needed for smaller R in order to keep the spot 
velocity (and thus the Doppler factor) at approximately constant 
level. 

The estimated radius of the spot (not a fitting parameter) cor- 
rected for orientation and gravitational effects varies between 3.0 
and 3.7 km for M = l.OM© and 1.6Mq, respectively. One can 
get some estimates of the size of the Comptonizing region from a 
rather large difference between the best-fitting Tcs and r obtained 
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Table 1. Best-fitting parameters 



M 


R 


R/rg 


i 


e 


r 


Tcs or a 


X^/dof 


(M0) 


(km) 




(deg) 


(deg) 












Model 1: 


Electron scattering slab, P'^ 


(p) OC /oo (tcs , m) 






1.0 
1.2 
1.4 
1.6 


c- , c+0.25 +0.45 
''■^'^-0.15 -0.35 
7 0+0.5 +1.1 
'•'-'-0.4 -0.6 
o Q+2.9 +4.5 
°-°-0.4 -0.7 
,o .+0.8 +1.6 
^•^■^-1.0 -1.7 


, 7r,+0.08 +0.15 
^•'^-0.05 -0.12 
1 04+0.14 +0.31 
^•"^-0.11 -0.17 

r, 1+0.7 +1.1 

^•^-0.1 -0.17 
n O+0.17 +0.33 
^•"-0.21 -0.35 


74+8 +8 
'^-20 -24 
09+0 +0 
^^-5 -18 

82^0 1° 
09+0 +0 

_2 


9r,+6 +10 
^'J-2.5 -3.5 
in + 1.5 +3 
^•^-1 -1.5 
„+0.5 +1 
^^-3 -4 
7+1 +1.5 
'-0.5 -1 


n ns+0'07 +0.09 
U.UO_Q Q4 _o.05 

n nfi+0-04 +0.08 

^'•^U-0.02 -0.03 
n ns+0-03 +0.05 

U.UO_g Qg _0.07 

n f.c:+0.02 +0.03 

U.UO_Q QJ _0.01 


r, a/l+0.22 +0.28 
'^•"^-0.16 -0.19 
n Kn+0.09 +0.24 
'^•'-"'-0.06 -0.10 
n tr9+0.05 +0.12 
'^•'^^-0.21 -0.22 
r, qn+0.04 +0.07 
'^•''■^-0.03 -0.03 


29.7/26 
35.3/26 
47.3/26 
61.4/26 








Model 2: 


P'^{^l) OC 1 - 








1.0 
1.2 
1.4 
1.6 


c- n+0.2 +0.4 
^■'J-O.l -0.2 
R 4+0.3 +0.8 
°-*-0.2 -0.3 

o q+0.6 +1.6 

''•■^-0.4 -0.7 
in s+0.6 +2.6 
l'-'-O-0.8 -0.9 


, f.7+0.07 +0.13 
^^•"'-0.03 -0.07 
1 70+0. 08 +0.22 
^•'°-0.05 -0.08 
1 Q7+0.15 +0.39 
l-y'-0.09 -0.16 
r, 9t; + 0.13 +0.47 
^•^0-0.17 -0.27 


09+0 +0 
°^-25 -32 
09+0 +0 
»^-8 -17 
09+0 +0 
»^-5 -13 

82^0 +Jr 


iQ+e +11 

^='-1 -2 
14+1.5 +3 
J^^-0.5 -1 
11+1 +1.5 
^^-1 -1.5 
c + l.O +1.5 
».0_0.5 -1 


n no+O-OT +0.09 
u-u»_o.oi -0.01 

n 1 .+0.03 +0.06 
'-'•^*-0.02 -0.03 
n 1 .+0.03 +0.04 
'-'•^4-0.01 -0.02 
p, 1C-+0.02 +0.04 
'^•^''-0.01 -0.01 


n fiQ+010 +0.12 

^^•"^-0.02 -0.03 

n 7.+0.O6 +0.10 

'^•'*-0.02 -0.04 

n sn+O-04 +0.08 

'-'•°'-'-0.04 -0.07 
n SQ+0-05 +0.08 
'^•"y-O.Ol -0.01 


28.1/26 
32.2/26 
38.7/26 
44.8/26 



in model 1. The values of Tcs are consistent with the optical depths 
obtained from the fits of the time-averaged spectrum with a Comp- 
tonization model for a slab geometry. For example, model COMPPS 
(Poutanen & Svensson 1996),^ gives Tcs = 0.5 for the electron 
temperature kT^ — 90 keV (which is not well constrained be- 
cause of the low signal above 100 keV). This leads us to a con- 
clusion that a large fraction of the black body photons reaches the 
observer without interaction with the hot Comptonizing medium 
(contradicting assumed homogeneous slab geometry). In order to 
get the "effective" r ~ 0.1 (see eq. 1141 ), about 40 per cent of the 
black body area should be covered by the Comptonizing medium 
with Tos ~ 0.5, while the remaining 60 per cent should emit as a 
pure black body (i.e. with no exp(— r//i) reduction due to scatter- 
ing). Thus, we conclude that the radius of the Comptonizing region 
^spot.c is about 2 km. 

The angular distributions of the scattered radiation (see dashed 
curve in Fig.|5ji) determined by parameters Tcb and a is also very 
similar for models 1 and 2. Model 2 is somewhat less restrictive and 
therefore gives generally better fits. The required emission pattern 
is certainly very much different from that given by the Lambert law 
and more radiation is escaping at intermediate angles a' ~ 50° 
than along the normal. It is this pattern which is responsible for the 
origin of a strong first harmonic (i.e. at double spin frequency) in 
the pulsar light curve at high energies (modified also by Doppler 
effects, see Fig.|4|and §|33J- 

In order to check the dependence of the results on the as- 
sumed geometry, we repeated the fitting procedure considering an 
extreme case with a point-like Comptonizing source at the cen- 
tre of the black body spot. The best-fitting parameters are simi- 
lar to the situation where both emitting regions coincide. This is 
expected since the size of the emission region should not play 
a large role since rspot = Rp < R- However, when estimat- 
ing the size of the spot, we assumed the soft emission is a black 
body, i.e. the observed color temperature is equal to the effective 
temperature. For the hydrogen and helium atmospheres of weakly 
magnetized neutron stars one expects that the color correction is 
/col = Tcoi/Toff = 1.6-1.8 (Zavlin, Pavlov & Shibanov 1996). 
The shift of the peak of the emission towards higher energies from 
the corresponding black body emission results from the fact that 
the outer layers of the atmosphere are cooler than the deeper layers. 

1 COMPPS is available at |ttp://ttp.astro.su.se/pub/juri/XSPEaCOMPPS| 
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Figure 5. (a) Observed normalized pulse profiles in the energy bands 3- 
4 keV (circles) and 12-18 keV (squares) and the best-fitting inodel light 
curves (solid and dashed curves) as a function of pulsar phase (two periods 
are plotted, zero phase is chosen when a spot is closest to the observer). 
A neutron star with the inass M = 1.4M0 and the radius R = 2rg = 
8.4 km and inclination i = 80° are assuined. The best-fitting parameters 
(model 2) are e = 11°, r = 0.16, a = -0.78 and x^/dof= 40.1/28. 
(b) The angular distribution of the intrinsic black body (solid curve) and 
Comptonized (dashed curve) fluxes /i/(/i) (here fi = coso') in the co- 
rotating frame normalized as J fil{fi)dii = 1. Only the range of angles 
between the dotted lines is actually observed, (c) The observed (crosses) 
and the model (solid curve) phase lags at the pulsar frequency relative to 
the 3^ keV band. 

Since the bound-free and free-free cross-sections rapidly decrease 
with increasing photon energy, one sees deeper and hotter layers 
at high energies. However, when the accretion takes place and the 
radiative shock forms (its existence is supported by the presence of 
the hard X-ray tail), the soft radiation can be produced by repro- 
cessing the hard photons from the shock. The temperature gradient 
of the atmosphere then can be opposite to that of standard neutron 
star atmospheres. Now the hot layers are at the top, softer photons 
are coming from the hotter region (since absorption cross-section is 
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large for small energies), while hard photons are coming from the 
cooler region below and the color correction can be smaller than 
unity. The actual vertical temperature dependence is however un- 
known. 

To investigate the influence of different possible color correc- 
tions on our best-fitting parameters, we consider two extreme cases 
of /col = 1-5 and 0.7 and fitted the data with model 2. The area of 
the "black body" spot now changes by a factor f^^i . For /coi = 1.5, 
the spot centre colatitude 9 increases by about 20 per cent com- 
paring to the case of no color correction. Since now the spot is 
larger, 9 has to increase to keep the same oscillation amplitude (see 
eq. 1221 ). The best-fitting stellar radius increases significantly only 
for the smallest considered mass of M = IM0 reaching 5.5 km, 
for M — I.2M0 the radius is 6.7 km, while for higher masses 
the change is negligible. In the case of fcoi = 0.7 the spot becomes 
smaller and also slightly decreases. At AI = IMq the best-fitting 
radius is now 4.8 km, while all other parameters remain the same 
within the errors. There is no change in the best-fitting parameters 
for higher stellar masses. This results from the fact that the light 
curve depends weakly on the spot size if it small comparing to the 
stellar radius (see eqs I21II22I ). 

One of the important assumptions that can affect our results 
is that the two emission components are assumed to be co-spatial 
or are at least co-centred. In our interpretation the Comptonized 
emission peaks earlier than the black body due to a different angu- 
lar distribution more affected by the Doppler effect. However, an 
additional phase shift is possible if the Comptonized emission re- 
gion is physically shifted on the stellar surface. We fitted the data 
with model 2 and the phase shift A(j) of the Comptonized emis- 
sion as an additional free parameter. A major effect of adding a 
new parameter is that the confidence interval for the stellar radius 
increases. Thus, for example, for M = 1.4Mq, the best fit with 
X^/dof= 34/25 was obtained for A(fi — 0.1 rad (a positive shift 
corresponds to an earlier arrival time) and other parameters very 
similar to those presented in Table Q while the constraints on the 
radius weaken: R = S.SIq'^ (90 per cent confidence limits). An 
_F-test gives F = 3.5 which corresponds to a 7 per cent chance 
that the decrease in was random (i.e. the improvement is not 
highly statistically significant). 

We show one of the fits to the pulse profiles with the fixed 
M, i, and R in Fig. |5^. The corresponding intrinsic angular dis- 
tribution of the radiation flux in the two spectral components is 
shown in Fig. [sJj. One sees a dramatic difference between the de- 
pendences of the black body and Comptonized components. The 
observed energy dependence of the pulse profiles and the phase 
lags at the pulsar frequency between different energies (Fig.|5j;) are 
reproduced reasonably well. A natural consequence of the model is 
that the phase lags change significantly at smaller energies together 
with the black body contribution to the total flux, while they reach 
a constant value at about 6-8 keV where the black body flux be- 
comes negligible (see also Gierlihski et al. 2002). Since in reality 
the angular dependence of Comptonization radiation is a function 
of energy, one would expect a weaker dependence of lags on energy 
even above 10 keV. 



5 DISCUSSION 

5.1 Constraints from the oscillation amplitude and the 
Doppler factor 

Let us consider first for simplicity the black body emission from a 
single spot. In that case, the Doppler effect does not change signif- 



icantly the light curve (compare dashed and solid curve in Fig.|4^) 
and thus it is mainly characterized by the amplitude A that can be 
estimated for a slowly rotating star using equation <19t . For a given 
R/vg, the amplitude in its turn depends on i and 6 only. Thus, at 
the i — 6 plane we can determine a curve that would correspond to a 
given amplitude A (see solid curves in Fig.|6^). One notes that this 
curve is symmetric around the i = 9 line since one can exchange 
i and 9 in the expression for A (see eqs I19H20I ) without affect- 
ing the result. If the black body emission is more beamed towards 
the spot normal because of a non-zero r, this curve will be some- 
what shifted down and left since A then increases. For example, the 
black body variability amplitude observed in SAX J1808.4— 3658 
is ^ 25 per cent (Gierlihski et al. 2002) and, correcting for r, we 
can determine the region at the i — 9 plane where it can be repro- 
duced for a given R/vg. 

The situation is different when we consider a less beamed 
emission pattern I (a) oc 1 + a cos a with a < (i.e. in our 
case radiation above ~ 8 keV). The oscillation amplitude is then 
strongly affected by the Doppler effect (compare dashed and solid 
curve m Fig.Et). The maximum of the Doppler factor and the spot 
projected area are shifted in phase by ~ 0.25 of the period and the 
peak of the emission is then shifted towards the phase where 5 has 
the maximum. However, the amplitude and the shape of the pro- 
file carry the information needed to determine a and the maximum 
Doppler factor (5max. Again at the i — 9 plane we can determine 
the curve of constant 5max (see dashed curves in Fig.|6t). Note, 
that since the Doppler factor (to the first order in (3) depends on the 
product sin i sin 9 (see eqs (9}ll II ). but not on i and 9 individually, 
this curve is almost symmetric around the i = 9 line. 

A very important point here is that the curves of constant A 
and 5max do not coincide (see Fig.|6|i). At the i — 9 plane thus there 
are at most two points (for a given star mass AI and radius i?) with 
the given values of these two parameters: points of the crossing of 
the curves of constant A and 5max. This leaves some ambiguity 
regarding the exchange of i and 9. It is possible, however, that the 
number of crossing points is one or even zero which would then 
mean that there is no physical solution with a given Doppler factor 
and the variability amplitude within the considered model. 

Let us now apply our simple arguments to the data on SAX 
J1808.4— 3658. The observed light curve implies a maximum 
Doppler factor about 1.015 and the variability amplitude A for a 
black body spot given by equation <19> of ^ 0.18 (actually ob- 
served higher value results from a sharper emission pattern with 
r ^ 0.1, see eq. 1141 ). The contours of constant A do not de- 
pend on the stellar mass but only on the radius measured in Vg (see 
Fig.lSJ)). We can now find a point at the i — 9 plane that satisfy 
both constraints. The curves of the constant A and 5max intercept 
in two points: at a large inclination and small 9 as well as at a small 
inclination i < 20° and large 9. 

The second solution is much less probable. First, we do not 
see any emission from the secondary antipodal spot. This implies 
that either it is blocked by the star itself (which in turn requires 
the stellar radius to be > 4rg for i ~ 15° and 9 ~ 60°), or it is 
blocked by the accretion disc. In the later case we have additional 
constraints. At inclinations smaller than ~ 50° an upper limit on 
the inner disc radius (i.e. the magnetospheric radius) required to 
block the secondary spot is _Rin < 27? (see lower dotted curve in 
Fig.|6|i). With such a small inner radius, the expected amplitude of 
Compton reflection component should be probably larger than the 
observed Q/2tt ~ 0.1 (unless the inner disc is highly ionised) and 
the shock should appear at uncomfortably large angle from the ro- 
tational pole if the magnetic field is a central dipole. On the other 
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Figure 6. (a) Contours of the constant variability amplitude A given by equation 1191 (solid curves) and the contours of constant maximum Doppler factor 
<5max (dashed curves) for a black body spot (v = 401 Hz, M = 1.4M0 and R = 2rg = 8.4 km). The dotted curves define the i — 9 relation where the 
photon trajectories from the antipodal spot (that without the disc would reach the observer) cross the disc at the minimum distance of 2R and the corotation 
radius Rco- These curves are computed using approximate photon trajectories from Beloborodov (2002). The inner disc radius smaller than 2R (or Rco) 
would block the antipodal spot from the observer in the area above the corresponding curve. 

(b) Contours of the constant variability amplitude A = 0.18 (solid curves) computed using equation 1191 for different stellar radii measured in rg units and 
the interception points of these curves with the contours of constant iSmax = 1.015 for different stellar masses: A/ = I.OMq (squares), 1.2Mq (triangles), 
1.4M0 (circles), and I.6M0 (diamonds). 



hand, at inclinations i > 70° the antipodal spot is blocked by the 
disc if the magnetospheric radius is smaller than the corotation ra- 
dius Rco = 31(M/1.4Mo)^''='(P/2.5ms)^''^ km which is also 
required for accretion to take place. Second, small i is inconsistent 
with the X-ray (Chakrabarty & Morgan 1998) and optical (Giles, 
Hill & Greenhill 1999; Homer et al. 2001) modulations at the bi- 
nary period, and finally, it is also inconsistent with the lower limit 
i > 28° obtained by Wang et al. (2001) based on the modelling of 
the optical/IR emission with the X-ray irradiated disc. We conclude 
that large inclinations i > 50° and small 9 ~ 5°-20° are preferred 
to small i and large 6. 

For small stellar masses M — I.OMq, the interception point 
for given A and 5max fixes the possible inclination i between 55° 
and 80° depending on the assumed radius. The minimum possible 
radius still satisfying constraints on the inclination angle i < 82° 
is R = 1.8rg. For higher masses, the minimum possible radius (in 
Tg units) increases, reaching R = 2.5rg for M = 1.4M0, and 
for M = I.6M0 there is no solution whatsoever. Thus it is quite 
natural that the is worse for the high mass star (see Table0. This 
simple analysis also explains why the higher is the stellar mass, the 
higher the inclination should be in order to satisfy constraints on 
both A and (5max. 

5.2 Stellar masses and radii 

Our results (see TableQand Fig.|7) put constraints on mass of the 
compact object. For M = IM0 the obtained stellar radius of i? ~ 
5 km is too small to be consistent with any existing equations of 
state for neutron or strange stars. For masses higher than I.6M0, 
the fits are much worse than for smaller masses. Thus we conclude 
that the stellar mass is bound between 1.2 and 1.6 solar masses. 

We can also put constraints on the equations of state of com- 
pact stars (which of course could be much stronger if the mass of 
compact object were known). Thus, if M — I.2M0 the only pos- 
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Figure 7. Constraints on the radius of the compact star depending on the 
assumed mass are shown by circles with double error bars (for 90 and 99 
per cent confidence intervals) for the fits with model 2. For comparison we 
show also the relations between the compact star mass and its radius for 
several equations of state of neutron stars and strange stars. The labels for 
neutron star equations of state are BBBl and BBB2 - Baldo, Bombaci & 
Burgio (1997), BPAL12 - Prakash et al. (1997), SBD - Sahu, Basu & Datta 
(1993), while for the strange stars SSI and SS2 correspond to the MIT bag 
model (see Gondek-Rosinska, Kluzniak & Stergioulas 2003 for details) and 
DEYl and DEY2 to the equation of state of Dey et al. (1998) computed as 
in Gondek-Rosihska et al. (2000). All realistic equations of state He to the 
right of the causality line. The line R = 2rg is shown to guide the eye. 
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sible solution is a strange star (Chakrabarty 1991; Li et al. 1999; 
Ray et al. 2000; Gondek-Rosinska et al. 2000) with the equation 
of state derived by Dey et al. (1998). For the mass about 1.4M0, 
almost all known neutron star equations of state predict larger radii 
than our best-fitting value. The exception are the models similar to 
BPAL12 from Prakash et al. (1997) and the models with a kaon 
condensate (Glendenning & Schaffner-Bielich 1999; see Lattimer 
& Prakash 2001 for more details). However, the maximum mass for 
these models is ~ 1.5Mq which is smaller than the best estimate 
of the neutron star mass in Vela X-1 of 1.86 ± O.IBM© (Barziv et 
al. 2001). Thus we are inclined to conclude that our constraints to- 
gether with the minimum neutron star mass in Vela X-1 put a lower 
limit on to the neutron star mass in SAX J1808.4— 3658 of about 
1.5M0. Strange stars (MIT bag models) with maximum masses 
in the range of 1.5-1.7M0 would also fit into our constraints. At 
larger mass of Af = I.6M0 a strange star is still a possible solu- 
tion as well as a neutron star described by equations of state similar 
to that of Baldo et al. (1997), while the most stiff equations of state, 
such as for example that by Sahu et al. (1993), can be ruled out (see 
Fig.|7}. 

Assuming alignment of the orbital angular momentum and 
that of the compact object, we can estimate the mass of the com- 
panion star Mc from our lower limits on the inclination of the 
rotational axis (see Table 0. Using the measured mass function 
= (Mcsin if/(M + Mcf = 3.78 x 1O"'^M0 (Chakrabarty 
& Morgan 1998) we find that Mc = 0.041 ± 0.002, 0.0445 ± 
0.0015,0.0478 ± O.OOO5M0 for M = 1.2, 1.4, 1.6M0, re- 
spectively. The companion radius is about 0.1 Si?© (Bildsten & 
Chakrabarty 2001). 

The mass measurements of radio pulsars in binary systems 
and their neutron star companions revealed their very narrow dis- 
tribution with the mean M ~ 1.35M0 and a very small dispersion 
(Thorsett & Chaki^abarty 1999). No evidence for extensive mass 
accretion was found. In order to spin up a star with the moment 
of inertia ~ 10*^ g cm^ to the frequency of 400 Hz one needs 
to accrete just about AM ~ O.IM©. Thus if the initial mass is 
within the same distribution one does not expect the mass of the 
compact object to exceed 1.5-1.6M0. For a mass of I.4M0, our 
best-fitting radius (with model 2) is i? = 8.31q'4 km (90 per cent 
confidence interval). Such a radius is consistent with the upper 
limit R < 7.5(M/M0)^/^ km obtained from the fact that pul- 
sations are observed at fluxes differing by a factor of 100 (Burderi 
& King 1998; Li et al. 1999). Here one assumes that the magnetic 
and the ram pressure from the accreting material are equal at the 
magnetospheric radius which is larger than the stellar radius R and 
smaller than the corotation radius Rco- For M = I.5M0, our best 
fit J? = 9.5 ± 0.5 km is already inconsistent with the given up- 
per limit. This constraint, however, is more relaxed if the magnetic 
field does not have a dipole structure and/or the accretion disc is 
not gas-pressure dominated (Psaltis & Chakrabarty 1999). 

5.3 Emission mechanism 

The angular dependence of the hard radiation that is needed to pro- 
duce the observed pulse profiles above 8 keV is consistent with that 
expected from the electron scattering dominated slab with Thom- 
son optical depth t^s ~ 0.3-1. It would be natural to assume that 
the emission originates in a plane-parallel shock at the surface of 
the compact star (Basko & Sunyaev 1976). The observed spectral 
slope r ~ 2 (the energy spectrum E~^) is expected from bulk mo- 
tion Comptonization in a radiation dominated strong shock (Bland- 
ford & Payne 1981). Most of the papers on this subject neglect gen- 



eration of new soft photons by reprocessing the hard radiation at the 
neutron star surface which would significantly soften the emergent 
spectra. Thus, in order to reproduce F ~ 2 slope, thermal Comp- 
tonization should play an important role (Lyubarskii & Sunyaev 
1982). As discussed in Gilfanov et al. (1998), there should exist a 
mechanism that adjust the spectrum so that it does not change much 
when the accretion rate varies by a factor of 100. Reprocessing of 
the hard photons produced in a shocked region into soft ones that 
cool the shock could be such a mechanism. This two-phase struc- 
ture in the energy balance is known to stabilize the spectrum at 
the observed values when the optical depth is of the order of unity 
(Haardt & Maraschi 1993; Stem et al. 1995). The electron temper- 
ature adjusts itself to the variations of the optical depth (and lumi- 
nosity) so that the product of the electron temperature and optical 
depth is approximately constant. 

Variation of the accretion rate that is the main candidate for the 
aperiodic variability observed in SAX J1808.4— 3658 (Wijnands 
& van der Klis 1998b) thus can cause the changes in the optical 
depth, but not in the spectral shape. However, the relative amplitude 
of the black body and Comptonized components can vary. Indeed, 
at small luminosities the black body seems to the more prominent 
(Gilfanov et al. 1998) which would correspond to a smaller optical 
depth of the Comptonizing layer. A larger relative normalization of 
the black body component should result also in a higher variability 
amplitude which in fact was increasing from 4 to 7 per cent (rms 
amplitude) over the duration of the 1998 April outburst (Cui et al. 
1998). 

There could be additional effects related to the electron- 
positron pair production (by photon-photon interactions) which 
becomes important when the compactness parameter I = 
Lar /{repot.c nieC^) is larger than unity. For the luminosity of 
L ~ 3 X 10'^'' erg s~^ (which is a characteristic luminosity dur- 
ing the April 1998 outburst of SAX J1808.4-3658, see Gierlinski 
et al. 2002) and the emission region size of rspot.c ~ 2 km, the 
compactness is about 400 and it could be much larger if the region 
has a smaller height. Thus, the data are consistent with the pres- 
ence of pairs which could be responsible for producing the optical 
depth in the range t^s ~ 0.3 — 1 depending on the compactness 
and for fixing the electron temperature in the range 50 — 100 keV 
working as a thermostat (Stem et al. 1995; Poutanen & Svensson 
1996; Malzac, Beloborodov & Poutanen 2001). If the hai'd radia- 
tion is slightly beamed towards the neutron star surface due to the 
bulk motion, the amount of soft photons produced would be larger 
than from an isotropic source and there will be no problem in re- 
producing the observed spectrum that satisfies the energy balance 
(cf. Gierlinski et al. 2002). 

The radiation escaping from the shock can be affected by cold 
plasma in the accretion column. Assuming the area of the emission 
region of iir ^pot.c ~ 12 km^ and the accretion rate M = 2 x 10^® 
g s^^ (needed to produce the observed luminosity at efficiency of 
0.15) one can estimate from the continuity equation the charac- 
teristic optical depth of the plasma over the distance equal to the 
spot radius r = n^aTrspot.c ~ 0.8. (At higher distances from the 
star, the magnetic field lines diverge and photons can escape freely.) 
This means that the flow can affect the radiation escaping along the 
normal to the spot and the radiation can be blocked from the ob- 
server at certain phases. Since the area emitting the black body ra- 
diation is larger than the cross-section of the accretion column, one 
would expect that the effect of shadowing is smaller for the black 
body. Possibly, this effect is observed in the light curve where the 
hard radiation has a break at (j> — (see Fig.|5^). We may also see 
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this effect from the best-fitting angular distribution (see Fig. (Sj)) 
which shows a drop along the normal at p = cos a' ~ 1. 

5.4 Origin of soft lags 

Our results confirm the proposal by Gierlinski et al. (2002) that 
the main physical reason behind the soft lags is a two-component 
nature of the spectrum. A natural consequence of this model is a 
break in the phase lag at the energy where the black body compo- 
nent vanishes, i.e. around 8 keV (see Fig.|5j;). We showed here that 
the black body and the Comptonized components having different 
angular distribution generate light curves which are affected by the 
Doppler effect in a different way. The Doppler shift is an important 
ingredient but not the major cause of the soft lags. In our model, 
the Doppler effect comes into play only through its influence on 
the variability amplitude of the bolometric flux. 

On the other hand, Ford (2000) and Weinberg et al. (2001) pro- 
posed a model where the lags in SAX J1808.4— 3658 are produced 
by the Doppler effect alone through its influence on the monochro- 
matic flux. In both papers a black body emission from the spot is as- 
sumed. When the spot is moving towards the observer, the Doppler 
effect shift the spectrum towards higher energies, while motions 
away from the observer softens the spectrum and generate a deficit 
at high energies. Thus, higher energy photons, in the tail of the 
black body distribution, arrive only at phases when motion is pre- 
dominantly towards the observer, while the flux of photons at the 
peak of the black body distribution has a maximum together with 
the projected area, i.e. quarter of the period later. This results in 
hard leads (or soft lags) as well. As discussed in Gierliiiski et al. 
(2002), there are a number of problems with this proposal. First, 
the photons in the black body tail contribute very little to the to- 
tal observed spectrum above 8 keV which is dominated by the hard 
Comptonized component. Thus, any soft lags in the black body flux 
would be completely invisible in the total flux. Second, a break in 
the phase lag energy dependence does not have any physical expla- 
nation in that model. 

Even stronger Doppler boosting is expected when photons 
from the spot are scattered off a much more rapidly rotating in- 
ner disc (Sazonov & Sunyaev 2001). The scattered radiation will be 
leading or lagging the incident spot radiation depending on whether 
the disc is co-rotating or counter-rotating. For a black body radi- 
ation from the spot, this would translate into soft lags or leads, 
respectively. However, since the spectrum from the spot is closer 
to a power law, no phase lags should appear (Chen & Shaham 
1989). A different angular distribution of the black body and Comp- 
tonized components as well as a non-coherent nature of scattering 
can affect this result. It could be an interesting problem for future 
studies. We believe that this model is not directly related to SAX 
J1808.4— 3658, since the spot is probably situated close to the rota- 
tional pole and the fraction of the spot photons scattered in the disc 
is negligibly small. 

Cui et al. (1998) interpreted soft lags as due to Compton 
down-scattering of intrinsic hard photons with energies above 10 
keV in a cold electron cloud of Thomson optical depth Tcb ~ 10. 
The softer photons (scattered more times) escape from the medium 
later producing soft lags. As the authors pointed out, the hard pho- 
tons should escape to the observer directly through some kind of a 
hole in the cloud since otherwise a spectral cutoff should appear at 
~ 511/res ~ 5 keV (e.g. Sunyaev & Titarchuk 1980; Lightman, 
Lamb & Rybicki 1981). This model, however, predicts a signifi- 
cantly reduced variability in the softer energy band (Chang Sl Ky- 
lafis 1983; Brainerd & Lamb 1987) contrary to what is observed. 



6 SUMMARY 

The two main spectral components, black body and Comptonized, 
are clearly identified in the time-averaged spectrum of SAX 
Jl 808.4— 3658. The pulse profiles show significant energy depen- 
dence and the softer photons lag the harder ones. Such a behaviour 
can be reproduced if the hard Comptonized emission peaks earlier 
than the black body emission (Gierlinski et al. 2002). The black 
body component shows about 25 per cent (peak-to-peak) variabil- 
ity amplitude and the Doppler effect (producing a ~ 7 per cent 
variability) does not change its pulse profile much (Fig.|4}. The 
Comptonized emission from an optically thin slab with a broader 
("fan"~like) angular distribution has an intrinsically smaller vari- 
ability amplitude. The Doppler effect then shifts significantly the 
peak of the pulse producing soft lags. 

By modelling the observed pulse profiles, we obtained an up- 
per limit on the mass of the compact object M < 1.6Mq. We 
also put the lower limit of M > I.2M0 on the stellar mass since 
for smaller masses the obtained radii are too small to be consistent 
with any published equation of state for neutron or strange stars. 
This also constrains the inclination of the system to be larger than 
65°. For the masses I.4M0 < M < 1.6Mq, the best-fitting stel- 
lar radii R = 8-1 1 km are consistent with those given by the neu- 
tron star equations of states as well as some strange star models, 
while for A4 ~ 1.2Mq with that given by the equations of state for 
strange stars only. 
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